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We first examine tiie microstructure of a cosmic string endowed witii two simpie Abelian currents. 
Tliis microstructure depends on two state parameters. We tfien provide tlie macroscopic description 
of sucii a string and show that it depends on an additional Lorentz-invariant state parameter that 
relates the two currents. We find that in most of the parameter space, the two-current string is 
essentially equivalent to the single current-carrying string, i.e., only one field condenses onto the 
defect. In the regions where two currents are present, we find that as far as stability is concerned, 
one can approximate the dynamics with good accuracy using an analytic model based on either a 
logarithmic (on the electric side, i.e., for timelike currents) or a rational (on the magnetic side, i.e., 
for spacelike currents) worldsheet Lagrangian. 

PACS numbers: 98.80.Cq, 11.27.-(-d 



I. INTRODUCTION 

Among all possible topological defects [TJ [2 [5] that 
may have appeared as consequences of early phase 
transitions, only cosmic strings may still be compatible 
with cosmology. They generically arise [U [5] in the 
hybrid inflationary models [6j |7j that are implemented 
in Grand Unified Theories [8', [9l [10] and are an un- 
avoidable consequence of Grand Unified Theories and 
supersymmetry breaking when there exists a low-energy 
unbroken i?— parity. More recently, cosmic strings have 
been the subject of renewed interest because it has been 
realized that such configurations could also exist as 
cosmologically relevant solutions in the context of string 
theory where D3/D3 and D3/D7 brane inflationary 
models lead to the formation of D-strings at the end of 
the expanding phase. It has been suggested |TT] that 
these objects correspond to D-term cosmic strings in 
supergravity theory. 

In many instances, e.g., in a supersymmetric the- 
ories [I2J, the Lorentz invariance that is typical of 
Nambu-Goto strings [131 [T?] used in cosmological 
settings [13 US HTj HHj US UHl HD, can be bro- 
ken by the appearance of currents flowing along the 
worldsheet [22] |23]. This effect is known to halt 
cosmic string loop decay caused by dissipative effects 
and yield new equilibrium configurations named vor- 
tons [24, i25l [261 [23 EH [29] . The density of these vortons 
is tightly constrained by the value of the normalized 
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density fl = p/pcrit today and by primordial nucleosyn- 
thesis |30l I31j . A non current-carrying string network 
also has a characteristic energy scale rrics which is 
constrained directly by the time stability measurements 
of binary pulsars |32j and indirectly by the observation 
that the Cosmic Microwave Background (CMB) [33] is 
dominated by the signal produced by the amplification 
of primordial density fluctuations during inflation [34j . 

Vorton (equilibrium) states can be destabi- 
lized PSI 1361 i37_| through various mechanisms, e.g., 
with a coupling of the Higgs field to an electromagnetic 
current [M] [39] [40], with the existence of shocks or in 
the presence of high curvature regions [41j [12]. The 
vortons are described by means of various possible 
internal equations of state [131 [HI US] US] HZ] which 
stem from a detailed study of the internal microscopic 
structure of the vortices that arise from the coupling 
of the string-forming Higgs field to a bosonic current- 
carrier [IHl HI] [SU]. More recently, vorton instabilities 
have (to some extent) been confirmed in large field 
theory simulations [STj . 

The macroscopic formalism developed by 
Carter [13] [HI US] HB] and that describes general 
p— branes embedded in a n-dimensional spacetime can 
be used to study the dynamics of current-carrying 
cosmic strings. As it stands, the formalism applies to 
the case of a single current for which a single state 
parameter provides a complete description of the string. 
However, it cannot be used for a system described 
by two or more state parameters as is the case when 
the Higgs field couples to fermions [5^. It gets worse 
when massive modes are present in the spectrum [53], 
because the current induces qualitatively different 
macroscopic properties when the state parameters are 
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changed. There is no complete description of a general 
"many-current" worldsheet and only the case of a "cold" 
superconducting current coupled to a "hot" entropy 
current appropriate for the "warm" string model has 
been treated i54j (see however the objections concerning 
the use of the word "superconducting" string in this 
reference) . 



products of all pairs of distinct currents. This number is 
equal to 1 when N = 2. In Sections [V] and VII we also 
consider the internal stability of a string endowed with 
N condensates (A'' = 2 in Section VII I and show that, 



contrary to what the three-dimensional case suggests, 
no new instability is predicted when several currents are 
involved. 



The two-current case provides the one-dimensional 
analog of the ordinary 3D Landau superfluid model |55j . 
It has been shown that hydrodynamical systems with 
more than one current present new |56j instabilities, 
as, e.g., the two-stream plasma instability [571 EH] (see 
Ref. [59] for a pedagogical presentation of this issue) . In 
order to test the stability of a vorton state against per- 
turbations in more than a single current, we provide in 
this paper one more step towards a complete formalism 
for the many-current carrying cosmic string. For that 
purpose, we consider in Section [IT] a specific model with 
two scalar fields trapped in the vortex core, leading to 
two conserved currents, and thus three boost-invariant 
state parameters. Being strictly local, the solutions of 
the field equations in the microscopic theory for the 
Nielsen-Olesen vortex ansatz [HD] depend on only two out 
of the three state parameters, respectively proportional 
to the amplitude of the two independent currents. These 
two state parameters therefore completely determine the 
string structure at any point along the worldsheet. On 
the other hand, the integrated quantities introduced in 
Section |IV| and the worldsheet dynamics described in 
Sections |V| and [Vll| exhibit an additional dependence on 
a third state parameter, a quantity proportional to the 
scalar product of the two currents. In fact. Section |V] 
deals with the worldsheet description of a string endowed 
with a set of N condensates, in which case there exists, 
in addition to the usual N state parameters, an extra 
set of N{N — l)/2 state parameters given by the scalar 



Finally, it is known that single current-carrying strings 
can be described by a macroscopic Lagrangian depending 
on a single string state parameter [47, 49J. In the case at 
hand, and since the field equations only depend on two 
state parameters, one could think of similarly describing 
a two-current string by means of a sum of Lagrangians, 
one for each current. In Section VIII| we compare the 
dynamics of the string as given by the fully interacting 
theory of previous sections to the dynamics of this sum 
of individual Lagrangians that describe non-interacting 
fields and that provide the possibility to perform a fully 
analytic treatment of the string physics. In principle, the 
two approaches do not describe the exact same physics, 
even when the coupling term of the fields in fully inter- 
acting theory is set to zero, but our results exhibit sat- 
isfactory agreement between the two for zero and even 
for weak coupling offering the possibility to study strings 
with several weakly coupled field condensates in a fully 
analytic way. 



II. TWO COMPLEX SCALAR FIELDS MODEL 

In order to introduce a current-current coupling in a 
vorton state, we couple a Higgs field H, charged under a 
broken local U{1) symmetry with associated gauge field 
and charge q, to two uncharged complex scalar fields, 
<i> and S, with global U{1) symmetry. 



£ = -- {D,H)^ {D,H) - -C^.C'- - - {\H\^ - rff - -d^<^*d^^<^ ~ ^d^^*d^^ - V {H, E) . (1) 
The interaction potential for (f> and E is given by 

V (H,^, E) ^ + + \ - V') iUm' + + ^A,!*!" + ^A.|S|4 + ||<i>nEp. (2) 

I 



The kinetic term for H reads Dfj^H = (9^ -I- igC^,) iJ 
and = d^Cu ~ d^Cf^. The vacuum is defined as 
usual as the minimum of the potential, i.e., through 
5V/5H = 5V/5^ = 6V/5T. = 0. This leads to a system 
of three cubic equations for the field amplitudes \H\, |<I>| 
and |S|. Although these equations are soluble in general, 
we shall only be interested in the regions of parameter 
space for which the absolute minimum is at |(i?)oP = if' 



(breaking the 1/(1)'°'^'*' ymmetry) and $ = E = 0. Such 
a choice is made implicitely with the potential written 
as in Eq. (|2|. 

There exist two additional [/(l)siobai gynimetries, 
carried by <i> and S, that independently leave the 
vacuum unchanged. Before symmetry breaking, 
one therefore has [/(l)'"'^'^' x t/(l)gi°bai ^ ^(j^)giobai^ 
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In the case = fa- and = = g, one has 

[/(;L)global ^ [/(j^-jglobal _^ 5'{7(2) global^ 

At zero temperature, the theory stemming from Eq. ([T]) 
admits vortex-hke solutions. In order to study locally 
cylindrically symmetric cosmic string configurations, the 
reference frame is chosen locally aligned with the vortex, 
and the z axis is defined to lie along the string. We denote 
by r and 9 the usual cylindrical coordinates centered on 
the string location and express the solution using the 
Nielsen- Olesen [50] ansatz 

H{x") = /i(r)e"^ and = ^^""^ Sl, (3) 

for a string with winding number n, h(0) = 0, and 
(5(0) = n. Far from the vortex, one recovers the vacuum 
state defined by h{oo) — \{H)q\ ~ rj and Q{oo) — 0. 

In this solitonic background, the breaking of 
^(]^)giobai ^ fj(2^)giobai^ dynamically generated as the 
bosonic fields $ and E condense in the string, can lead to 
the appearance of currents along the vortex. To see this, 
we repeat the analysis, first presented in Refs. [231 US] 
and based on a perturbative expansion in the fields $ 
and S which are assumed small to begin with. The field 
equations derived from Eq. ([T]) can be written in the 
form of two-dimensional time-independent Schrodinger 
equations as follows. Given <I> = (j){r, 6')e*^'** and S = 
a{r, 0)e'^''* and neglecting nonlinear terms in the fields, 
one has 

~AcP + V^.p = El(t>, (4) 
-Aa + Va<J = Ela, (5) 

with 

V^^ f^{h'-7f) + ml, Uih^-r,^) + ml. (6) 

If these potentials are both negative definite, there exist 
bound states for both fields, i.e., solutions with < 
and < 0, leading to instabilities in these fields, and 
therefore to condensates. If the vacuum masses of $ and 
E vanish, this is certainly the case; it holds true as well 



in a neighbourhood of (m^,TO^) = (0,0) provided the 
constraints 

UV^ >ml, Uv^ > ml (7) 

are imposed. Although these conditions have been de- 
rived from a dynamical analysis, they can be recovered 
by examining the vacuum condition SV/S^ — SV/ST. — 
at the string location, i.e., setting H = 0. Assuming nei- 
ther $ nor S vanish, this gives 

A^|<i>P + .g|E|2 = Uv'-^l, 
<7|<i>|2 + A,|E|2 = fa^'-ml, 

which is only possible in the range of parameters defined 
by Eq. Q. As we now make clear, although these con- 
straints are necessary they are by no means sufficient. In- 
deed, when one field condenses, inclusion of its nonlinear 
potential term ((7/2)|<I>p|Sp in the action modifies the 
other field Schrodinger equation in a way that cannot be 
determined analytically. This correction is positive def- 
inite so that a strong coupling between current-carrying 
fields tends to drive one of the condensates to zero. In 
this work we focus our attention on the effect of this cou- 
pling on the string's energy density, tension and stability 
and we further restrict attention to the weak to moderate 
coupling case and impose the condition 

A^A, > (8) 

on the quartic coupling constants. 



III. TWO-CURRENT-CARRYING 
CONFIGURATIONS AND NUMERICAL 
SOLUTIONS 

For the cylindrically symmetric forms 

$(x") = (t,{r)e''^'>' = (f,{r)e'^'^*'-''*'\ (9) 
S(a;") = CT(r)e*'^'- = CT(r)e''("'<'*-'^'"^\ (10) 

together with ([3]), the field equations read 



d^h 1 dh 

dr^ r dr 

d^Q 1 dQ 

dr^ r dr 



d^(/) 
dr^ 

dV 

dr^ 



1 d(j) 
r dr 
Ida 
r dr 



2 

2 1,2 



h, 



[wcr + U{h^ - rj^) + ml + y^acr"^ + g(j)^]a, 



(11) 
(12) 
(13) 
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where the two real state parameters 
W4> = kl~ Lul and 



_ 7 2 2 



(14) 



are Lorentz scalars. Knowledge of these two parameters 
is sufficient to fully determine the microscopic structure 
of the defect. Since the currents are vectors, an additional 
Lorentz scalar can be built out of the vortex fields. It 
represents the relative position of the two currents and 
arises in the two integrated functions (energy per unit 
length and tension) that are necessary to describe the 
string dynamics completely. 



A. Dimensionless quantities 



The coupled equations to \1A\ can be solved using 
a Successively Over Relaxed (SOR) method, a procedure 
appropriate to this category of boundary condition non- 
linear systems |48l |49j |50] . This method is presented in 
full detail in Ref. pTl. Defining the dimensionless dis- 



tance variable 

p = VXrjr, 

and the rescaled field functions X, Y and Z through 



we renormalize the coupling constants as 

.,2 



~2 _ T 



and 



2 r 7 4 

A,r;2' '''-AA.77^' " AA,r;4 ' 



(15) 



(16) 



(17) 



the index i standing for either (jj or a. In terms of the 
dimensionless variables and couplings, the action reads 



S = - 



TTTj 



dxy x'Q' 1 

dp ) p2 q2p 



dg 
dp 



+ ^(^2-1)2+^^ 



dr 
d^ 



)(X2-l) + ;^70i^4+;^7.^^ + 



dZ 

dp 



pdp. 



(18) 



in which we have used the dimensionless form of the state 
parameters. 



Wi 



(19) 



with i = 0,cr. The constraints ([t]), when expressed 
in terms of the dimensionless parameters, simply read 

A > Ir- 

The coupling parameter space was explored by fixing 
ai, Pi and q and varying Ai and g for light-like currents, 
i.e., for Wi = 0. Three typical solutions are presented 
in Fig. [T] These solutions were obtained with the fol- 
lowing boundary conditions^: Y'{Q) = Z'{0) — Y{oo) — 
Z{oo) = 0, while the vortex itself satisfies X{0) = 0, 
X{oo) = 1, Q(0) = n = 1 (we restrict attention to unit 
winding numbers strings in what follows) and Q'{0) = 0. 



^ For now on, a prime will represent a differentiation with respect 
to the rescaled distance variable p. 



B. The state parameter space 

We now discuss the allowed range of values for Wi. We 
first consider the solution near the string core. Assuming 
two condensates to be present, we expand both functions 



Y ^yo + y2P 



Z = zn + Z2P^ + • • • , (20) 



Using ( 20 ) in the field equations yields 



W^<(3^-"^^. (21) 

From Eq. ([t]), we note that (making use of the one- 
current terminology [49] ) this is constraint on the 
"magnetic" side of the current, i.e., a constraint on the 
possible range of spacelike currents having w > 0. 

For p ^ oo, the field equations read 
1 



-Y' 



Z' 



(w0 + 70) Y, 
(wa + 7a) Z, (22) 



and their solutions are expressed in terms of Bessel func- 
tions. As in the one-current case, there are phase fre- 
quency thresholds above which the currents radiate away 
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from which the microscopic structure can be derived 
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FIG. 1: Dimensionless field profiles obtained using = 1.8 x 
10-^ = 1.4 X 10-^ = 1.0 X 10-^ = O.S X 10-^ 
70 = 5.5 X 10"*, 7a = 4.5 x 10"'', = 1.0 x 10-\ w^, = 
Wa — X — and finally from top to bottom ~ Qi ~ 0.0, 
<f = 92 = 0.05 707a and = gl = 0.5 7^,7^. 



must satisfy —mf < Wi < fifj 



These conditions 



are necessary to confine the current, but not sufficient 
to ensure there is indeed one. Indeed, because of the 
non-Hnear interactions between the fields, even when 
these conditions are met, there exists the possibility that 
the total energy be minimized by a vanishing current. 



IV. EQUATION OF STATE 

In this section, we determine the macroscopic, i.e., in- 
tegrated, quantities from which the string dynamics can 
be derived. These are the total currents associated with 
^{x") and T,{x°') and the eigenvalues of the stress-energy 
tensor. 



A. Currents 

A conserved current is associated with each [/(l)siobai 
symmetry. Given ^ and (10), one has 



where 



(24) 



(25) 



The unit vectors and are respectively timelike and 
spacelike and are defined by 



and 



(26) 



in the system of coordinates {t, z,x,y}. Given this con- 
vention, the two state parameters Eq. ( 14 1 are 



w. 



g^'^d^ipidiyipi (no sum on i = a,t 



(27) 



What matters in a macroscopic formalism such as that of 
Refs. |43l HH HS] |46] is the Lorentz square of the currents 
integrated in the transverse direction. 



C = \/ CuC , 



(28) 



from the vortex. This translates into a constraint on 
timelike currents, 

Wi > -7^. (23) 
Both constraints imply that the two state parameters 



where 



(29) 



We denote the integrated currents hy d^, i — (j),a. For a 
straight and static string configuration, we then have 
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W0 



(j^rdr 



2x 



dr 



a^rdr 



(30) 



which involves, in addition to and w^, a third state 
parameter, given by 



(31) 



We note that, as mentioned in the introduction, although 
it does not enter into the microscopic description, x 
appears in the macroscopic description of the vortex 
dynamics. 

In the macroscopic formalism, the integrals of the 
scalar fields over the transverse section are especially 
meaningful. They read 

JC^ = —2n j (f>^{r)rdr and /Co- = — 27r J a^{r)rdr, 

(32) 

and correspond to the partial derivatives of the macro- 
scopic Lagrangian, 



C{wi 



(33) 



with respect to the currents. Given that the solution of 
the field equations are determined with the 4-dimensional 
Lagrangian, variations with respect to the parameters Wi 
yield directly 



K-i = 2 



dwi 



(34) 



The integrated current C can of course be re-expressed 
using dimensionless variables. Defining 



we then have 



7? 



(35) 



C = 



~ la 



Z^pdp 



l<t>la 



2iJ^^ { I Y^pdp 



Z^pdp 



1/2 



(36) 



r 



where the third dimensionless parameter x reads 



(37) 



B. Stress-energy tensor 

The stress-energy tensor can be obtained from the La- 
grangian through 



5C 



Sg' 



g^'C. 



(38) 



It has two eigenvalues corresponding to a spacelike direc- 
tion and a timelike direction of the worldsheet. Denot- 
ing by and the normalized, respectively timelike 
and spacelike, eigenvectors of T^'', and setting r|^^^ — 



—u^^u'^ + v'^v^, the first fundamental tensor of (whose 
mixed form is the projector on) the string worldsheet, 
the eigenvalues, namely the energy per unit length U 
and tension T, are obtained through the expression for 
the integrated stress-energy tensor |331 HH 251 133] 



(39) 



Unlike the case of a string with a single current, T^'' 
obtained in Eq. (38 1 is not automatically diagonal. We 



therefore define (in the worldsheet coordinates t and z) 
~ab _ Diag {— -1-1}, the string metric tensor, and ex- 
press the two-dimensional part of the stress-energy tensor 
as the sum of a diagonal part T"^ = —Afj°'^ with 



and a mixed non-diagonal part 



T 



,ab I B C 



ND \cBr 



with 



and 



C — 2n {k^ujfi,4>^ + ka-UJa<j'^) rdr. 



(41) 



B = n [{kl + ojD 4>^ + {kl + Lol) a^] rdr, (42) 



The eigenvalues of the stress-energy tensor are obtained 
through 



det iT"^ + h~f^) = 



(44) 



and we find that the energy per unit length and tension 
read 



u = a + ^/b-^-c^, 

T = A- v/B2-C2. 



(45) 
(46) 



(43) The quantity 



W(j, I (j) ''dr 



Wo 



cr^rdr 



rdr 



(7^ rdr 



(47) 



is obviously positive definite, so that U and T are well- 
defined. Using (32), this expression can be re- written 
as 



B'^ - C'^ ^ { ^WijjK.^ - ^W„ICa 



+ x^lC^JCo. (48) 



Note that the Nambu-Goto equation of state is recovered 
when B"^ ^C^. 

Both U and T are shown in Fig.[2j in the {f^ , z'^) frame. 
For negative and small positive values of w^_o , there are 
two condensates in the string. For larger positive values 
of either or w^, one of the two condensates vanishes 
and the physics reduces to the usual case: the tension 
diverges to negative values for Wi while the en- 

ergy density diverges to positive values (see Figs [3] and 
|4]). When there are two condensates in the string, there 
exists an additional divergence when both and re- 
spectively tend to and —"fa- In this limit, the quan- 
tity A diverges because it contains an integral in r over 
a sum of the squares of <i> and E in the potential term 
and they are themselves divergent [50,- When x = 0, 



\/ B"^ — on the other hand is a slowly growing func- 
tion in this limit because it is equal to an integral in r 
over the difference of the squares of <i> and S. Note that 
in Figs. |3] and [4] the full lines represent one dimensional 
profiles of U and T along and Wa for different values 
of Wo and respectively. The dashed lines represent 
the approximate values of U and T obtained from the 
analytic model |47l 149] which we compare to the numer- 
ical results in section VIII The timelike and spacelike 



eigenvectors Ua and Va are defined respectively through 
the eigenvalue equations 



T^^Mfc = -Uf]''^Ub, with 



and 



nb^-rr^vb, with fj^%,vb^+i 



rf^UaUb = -1, 



(49) 



(50) 



The canonical form ( 39 1 is recovered provided the four 
dimensional eigenvectors defined through Eqs. (49 1 and 



(50 1 are chosen to be the solutions of Eqs. (49 1 and (50 1 



The eigenvector thus read 



B + v/S2 - C2 



-1/2 




(51) 
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B - x/B^ - C2 



-1/2 



' B - y/B^ ~ C2 ' 
C 






(52) 



Again, given that B^ — > 0, the quantities under the 
square roots in and are positive definite, a result 
which will be useful below. In addition, the combination 
B^ — depends only on the state parameters and on 



the JCi's and is thus a Lorentz scalar. Therefore Eq. (48) 
holds in whatever frame. One may therefore work in a 
frame in which C 



gives the value of B. Note that in order to see that this 
limit is indeed well-defined in Eqs (51) and (52), these 



two expressions need to be rewritten in a slightly different 
form from the one presented here. In the limit C ^ 0, 
one then finds that they are nothing but the unit vectors 
and of Eq. ( 26 1 . This frame is the generalization 



0, in which case Eq. (48) simply of that frame for which, in the single current case, the 
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FIG. 3: U (upper curves) and T (lower curves) one-dimensional profiles along for Wa = —2.25 x 10"'' ~ 0.0 and 

2.25 X 10"* ~ 7i/2 (from left to right) and for gi, §2, ffs (from top to bottom). The solid lines correspond to U and T 
obtained numerically while the dashed lines were obtained from the analytic model |47l 149) discussed in Section VIII Note the 
^-coordinate defined by v ; 



in order to emphasize the neighborhood of iZ; = 0. 



phase gradient of the current carrier depends either on 
time or on space, but not on both, i.e., the boosted frame 
in which either the frequency or the momentum of the 
trapped scalar field is removed. 



C. Available ranges of variation of the underlying 
parameters 

There exists a finite w^-range within which two 
currents can appear. In this range JCi ^ {i = 1,2) 
and X is well-defined. The coupling between $ and S 
acts (nonlinearly) as a positive mass term for both these 
fields and quite generically (for any given set of masses 
and coupling constants) will reverse the condensation of 
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FIG. 4: Same as Fig. [s] but along w„ with = -2.25 x 10"* ~ -7^/2, 0.0 and 2.25 x 10"* ~ 7^/2. 



either of the two fields outside the appropriate w^-range. 
The string then behaves like a one-current-carrying 
string. 

The range of variation of x can be constrained using 
the spacelike or timelike character of the two currents. 
Setting generic two-dimensional normalized spaceHke and 
timelike vectors 5^ and respectively, with 



Stiil) = ipsinh7-|- z^cosh7, 
Tfih) = ipCosh7 + z^sinh7, 



(53) 



where 7 is a constant, the currents c* can be chosen 



proportional to either S"^ or T^, depending on their 
nature. If the two currents are of a different kind. 



with, e.g. 



spacelike and c}, timelike, we can take 



ci = ^m^Sf^i^j^,) and 4 = \/-w^7)i(7o-); lead- 
ing to a: = —yi'—w^Wa- sinh(70 — 70-), which provides 
no further restriction on the range. On the other 

hand, if both currents are spacelike, we have c'/f "' — 
^yw^Sfjl'ycf,) and c)j — ^/w^Spt{'^a), and we then find 
that X — ^W(jjWc cosh(70 — 70-), and is therefore larger 
than ^w^Wij. Similarly, if both currents are timelike, we 

set c^*'' = ^-w<^r^(70) and c^^^ = -WaT^,{^„), lead- 
ing to the result that x = — ^w^w„ cosh(70 — 70-), which 
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is now smaller than —^w^Wa. In short, for two currents 
having the same character, the range of variation of x is 
restricted to 



\x\ > ^W^Wo 



(54) 



One of the currents, cjj*''' say, can be lightlike. In this 

case, it must read c^*'' = TO^(ei^ + e'z^), where we 
have fixed the arbitrary normalization to the mass 
of the corresponding current-carrier for reasons of 
dimensions, while = e'^ = 1. If cj^^ is spacelike, i.e.. 



/waS^{'^a), we obtain x = e'm^^yw^e 



while if it is timelike, with c/, — ^/—'WaT|^{'^cr), the 
third state parameter is then x = —em^y/—Wij,e~'^'^ '^^ . 
Finally, both currents can be lightlike, with 



„(*) 



'm^{e<ptfj, + e'.Zfj,) and c 



leading to a; = -m^m^{€^ea - e^e^). 



We shall see in detail in Section VII A| that the range 
of X is further restricted by imposing the stability of the 
string under transverse perturbations. 

It is important to realize at this stage that the macro- 
scopic dynamics of a string on which two condensates 
can appear is more complicated than its one condensate 
counterpart. In particular, the number of degrees of free- 
dom is itself a dynamical variable: as either state param- 
eters and evolve along the string or with time 
(they both are in principle functions of the string inter- 
nal coordinates [3TJ US |331 HH US US]), the correspond- 
ing fields may switch back and forth between condensat- 
ing to non-condensating situations, with the consequence 
that the number of state parameters may jump discon- 
tinuously, being equal to either one or three (we are as- 
suming that the underlying parameters are such that for 
Wij, = = 0, both condensates are present). This seems 
to forbid any kind of macroscopic treatment such as pro- 
posed in Ref. gSl HH US 05] . We shall see however that 
because the string remains essentially a one-dimensional 
object, its classical stability can be investigated provided 
a generalization of the usual framework is made. This is 
was we do in the following sections. 



V. ELASTIC STRING DYNAMICS 

The dynamics of a current-carrying elastic string de- 
pends on the string's internal degrees of freedom but 
not on its geometry. It has been extensively studied in 
SH HH 135]. Because its Lagrangian depends only 
on the two state parameters Wi, a string with two con- 
densates can be studied using the same formalism. The 
purpose of this section is to extend the existing formal- 
ism for the dynamics of a single-condensate elastic string 
to an A^-condensate string. 



A. Preliminary geometric definitions 

The string worldsheet is defined as the two-dimensional 
surface swept by the string during its time evolution. 
It has timelike and spacelike directions associated with 
timelike and spacelike internal coordinates and 
The metric on the string worldsheet reads 



(55) 



where the subscripts '\a"i ••■ denote derivation 

with respect to the inverse metric is h°'^ . 

The embedding of the two-dimensional metric in four di- 
mensional spacetime is defined through the string world- 
sheet's first fundamental tensor as 



77- 



(56) 



Its mixed form 77^ is identified with the tangential pro- 
jector on the string worldsheet. The orthogonal projector 
is then defined through 



(57) 



With these definitions, it is possible to embed all fields on 
the worldsheet in four dimensional spacetime, provided 
that the 4D covariant derivative taken along the string 
worldsheet Vq is replaced by a new longitudinal covariant 
derivative which projects out the meaningless varia- 
tions transverse to the string worldsheet, 



(58) 



We also introduce the anti-symmetric fundamental tensor 
through the relation 



(59) 



Note that e is defined up to an overall sign. The curvature 
tensor, 



(60) 



is tangent to the worldsheet in its first two covariant 
indices and v and orthogonal to it in its contravari- 
ant index p. Since the projector 77''^ defines the space 
tangent to the string worldsheet, K^^'' is symmetric 
in its first two indices. This is the Weingarten iden- 
tity 1331 Si SUSS] 



[^lv 



p — 



= 0. 



(61) 



Finally, using two orthonormal basis vectors and 
tangent to the worldsheet and chosen timelike and space- 
like respectively, the fundamental tensors take the simple 
form 



(62) 
(63) 
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B. The elastic string model 

As stated above, an elastic string is described by an 
effective Lagrangian which only depends on its internal 
degrees of freedom. In the case at hand, there are N in- 
ternal degrees of freedom given by the N fields tpi associ- 
ated with the N currents living on the string worldsheet. 
The Lagrangian 

^ = ^x^J) (64) 

may thus a priori depend on any of the scalars 



(65) 



i.e., the full set of state parameters that form an 
N X N symmetric matrix. The Lagrangian will thus 
be written in the remainder of the paper as a function 
of the symmetric matrix x = iXij) denoted by C{x)- 

For instance, in the case of the two currents presented 
in the previous sections, Eqs. (27 1 and (31) show that 
and X0cr 



Xi 



X(y<p 



whereas Eq. ( 33 1 immedi 



ately shows that the Lagrangian of this particular model 
only depends on the diagonal entries Wi of the matrix 

The equations of motion of the elastic string can 
be obtained directly by varying the Lagrangian with 
respect to the worldsheet coordinates x^{£^°') and to 
the internal fields ipi. However, it is easier and more 
useful to write the dynamical equations as conservation 
equations because the physically meaningful unknowns 
are the string conserved currents , not the internal 
fields ipi. 

There are 27V + 2 independent degrees of freedom: 
six for the first two currents which define the string 
worldsheet through its tangent space and therefore 
satisfy the Weingarten identity (61 1, and then two for 
each of the other N — 2 currents which live in this same 
tangent space. 2N + 2 equations are therefore needed. 

The N conserved currents inside the string 

where 5^ j is the Kronecker delta, can be expressed in 
terms of the phase gradients through 



with 



(67) 



(68) 



Note that Eq. ( 68 1 generalizes Eq. ( 34 1 . These currents 



are an obvious choice for a first set of N conservation 
equations. Their conservation equations read 



V^(c'^) = 0. 



(69) 



The stress energy momentum tensor T^'' defined by Eq. 
( 39 1 or equivalently by 



5C 



satisfies the conservation equation 



/:?7m--/c*^(v^V0(v.v,), (70) 



(71) 



and provides 4 additional conservation equations. The 
first two are given by the transverse part of this relation, 
namely 



±/V^T^" = 0, 



(72) 



defined as its projection orthogonal to the worldsheet, is 
associated with the geometry of the string worldsheet. 
The remaining two equations are given by the longitudi- 
nal part of Eq. ( 71 1 , 



^/.''V^T'^'^ = 0, 



(73) 



and is associated with the internal degrees of freedom of 
the string. Finally, there exists an irrotationality condi- 
tion on each of the N gradient fields V^Vii 



(74) 



The vector V^(e'^'') being purely orthogonal (see 
Eq. (63l), these irrotationality conditions can be turned 
into conservation equations through an integration by 
parts, 



where 



V^«) = 0, 



(75) 



(76) 



The previous considerations result in the existence of 
27V conserved currents and a conserved tensor for a 
total of 27V -|- 4 conservation equations. These are 
two more than needed. This is because, as will be 
shown next, there exists a redundancy between the 
two longitudinal stress energy tensor conservation 
equations and the current conservation equations. In 
the well known case of one current only, i.e., for 7V = 1, 
the two current conservation equations are exactly 
equivalent to the longitudinal part of the stress energy 
tensor. On the other hand, in the general case, the nat- 
ural choice is to keep the 27V conservation equations only. 



We now show that the conservation equations ( 69 1 and 



( 76 1 , or equivalently ( 74 1 , imply the longitudinal conser- 



vation of the stress energy tensor. Rewriting the stress 
energy tensor as 

Tf"" = £77'"' - c''^(£%dp = £77^" - (e^pOc"', (77) 

we have 

?/,,pV,T'^'^ = Vp£ - Ep^c^'^V.df . (78) 
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Furthermore, the gradient of Eq. (64 1 reads 



from the other currents di^, and the symmetric matrix 



so that Eq. ( 78 ) becomes 



(79) 



(80) 



Finally, c^^rj'^^P is tangential and antisymmetric in its fi 
and v indices and is thus proportional to the antisym- 
metric tangential tensor e'"', 



Eq. (78 1 therefore reduces to 

?7MpV.T^'' = (£p.c'")(V,0, 
and as expected implies 

Vt.p^uT^" = 0. 

C. Equations of state 



It is clear from the previous section that the equations 
of motion can be reduced to the 2A'' + 2 conservation 



(81) 
(82) 
(83) 



equations (69 1, (72 1 and (75 1. As stated before, the cur- 



rents are the physically most meaningful quantities 
and, from their definition (76 1, yield the matrix of state 
parameters 

X^J = -ird,^dj,. (84) 
They also appear explicitly in N of the conservation 



equations (75 1 



On the other hand, the iV + 2 other equations do not 
depend explicitly on the 's. However, since they deter- 
mine the matrix of state parameters through ( 84 1 , which 



in turn determines the Lagrangian and its derivatives, it 
is certainly possible to express the stress energy tensor 
T^'^ and the other conserved currents as functionals 
of the unknowns d^^. These functionals are therefore 
additional identities that form a set of equations of state. 
In general, they will be determined by the underlying 
field theory model (i.e., a generalization to N currents 
of what was done in the first four sections of this paper) , 
and will characterize the current-carrying string model 
under consideration. 



The iV^ equations of state expressing 



(85) 



are necessary and sufficient to determine the conserved 
currents 



X given by (84) 



As for the stress energy tensor, it is a tangential, and 
thus two-dimensional, symmetric tensor which is there- 
fore determined by three parameters. The stress energy 
tensor is usually expressed in the diagonal form ( 39 1 



where and are a basis of respectively timelike and 
spacelike orthonormal vectors, whereas the eigenvalues U 
and T can be respectively interpreted as the string en- 
ergy density and tension. Thus, the natural choice for the 
three parameters defining the stress energy tensor are U , 
T, and the position of the basis {u^jV") in the tangent 
space. The latter is akin to a hyperbolic angle -0 with a 
given direction such as that of one of the currents di^, 
di^ for example. Note that this only defines the vectors 
u'^ and v'^ up to a sign, but this sign does not change the 
stress energy tensor ( 39 1 . There will thus be three extra 
equations of state expressing these three parameters as a 
function of the matrix x- In conclusion, there are iV H- 3 



equations of state which will be chosen here as ( 84 1 , (851 
and 

C/=C/(x), r = T(x), = 0(X), (87) 
where the hyperbolic angle is such that 

dif, cx cosh(V')u^ + sinh(V')w^, (88) 
if dif^ is timelike, and 

dif, cx sinh(V')u'' + cosh(i/')w^, (89) 
if it is spacelike. 

D. Duality 

In the case of one current, there exists a duality which 
allows to exchange the roles of the currents and 
[ISl HU HH US]- The existence of 2N conserved currents 
instead of the expected N suggests that there similarly 
should be a duality between the two sets of A'^ conserved 
currents c*'^ and d^. This is a duahty under the global 
exchange of the two sets of currents, not just of two given 
currents c*'' and . More precisely, we are looking for a 
dual Lagrangian C, of the same general form as £ (but 
function of different scalar fields ■0), yielding the same 
equations of motion, but such that the two sets of N 
conserved currents are exchanged. Using tildes to distin- 
guish all quantities derived from the dual Lagrangian, it 
is sufficient to require that 



(90) 
(91) 
(92) 



Using the notations of Eqs. (64 1, (65 1 and (68 1 as well 



(86) as the definitions (66 1 and (76 1 with tildes for the dual 
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model, Eq. (91) becomes 
which yields in turn that 



X 



(93) 



(94) 



where the right hand side is to be understood as matrix 
multiplications. Then, substituting Eq. (93) into (90) 
yields 

V^i^, e^'^/C.jVW^ = t,jK.^''V ^^u, (95) 



which in turn implies 



/c = k; \ 



(96) 



where the right hand side is again to be understood as a 
matrix inverse. 



Finally, the last equation ( 92 1 can be expanded using 



Eq. ( 70 ) and the previous expressions for JC and V-^i as 



The identity 



(97) 



(98) 



which can be derived for instance by using a pair of 
orthonormal vectors and Eqs. (63), enables to simplify 



Eq. (97) to 



(99) 



where in the latter expression, the matrix trace Tr has 
been taken. 

Thus a model with Lagrangian C {xij ) given by 
Eqs. (94), (96) and (99) has, by construction, the same 



2N conserved currents and stress energy tensor, and thus 
the same equations of motions as the initial model de- 
scribed by the Lagrangian C. However, these three equa- 
tions are not completely independent since one must sat- 
isfy consistency equations derived from Eq. ( 68 1 , namely 



(100) 



This equation can straightforwardly be derived from 
Eqs. (Iggl) and (l94|. 



Thus, the model described by the Lagrangian C {x) 
derived from the initial model by Eqs ( 94 1 and ( 99 1 pro- 



duces the exact same elastic string dynamics. This is 
precisely the duality found previously SH 113 SB] for 
an elastic string with a single current. 



VI. THE ELASTIC DOMAIN 

The elastic A^— current-carrying string model intro- 
duced in preceding sections is obtained from the field 
theory solution of a straight and static string solution 
and further assuming that, at the macroscopic level, 
the string is locally straight and static in a rotated and 



boosted frame in which and defined in Eq. (39) 



are aligned with the t and z axes. This is usually an 
excellent approximation since the string curvature is 
of the order of the Hubble scale, while its thickness 
(to which the latter curvature must be compared to 
evaluate the "straightness" of the string) is of the order 
of the Compton length of the lightest particle involved 
(typically a current carrier). 

However, this approximation is only meaningful if the 
static straight solution is dynamically stable, because 
otherwise, the string curvature, however negligible to 
begin with, will drive the string away from the straight 
and static solution. The conditions under which the 
string remains straight and static will define the domain 
of elasticity of the string. Evaluating string stability 
at the field theory level is a complicated task. Instead, 
we shall determine whether the static straight string 
solution is dynamically stable in the elastic string model. 
If it is, the description is self-consistent. If it is not, 
a more detailed calculation must be performed at the 
microscopic level. 

In the straight and static solution, the currents c?f are 
constant along the string worldsheet. This implies that 
the other currents c'^ and the stress energy tensor T^"^ 
are also constant. This ensures that the equations of 
motion (69), (72) and (75) are satisfied. 



Here, we shall take the static timelike Killing vector 
and the static spacelike Killing vector z** introduced 
in ( 26 ) to describe the time coordinate and the direction 



along which the string lies. The stress-energy tensor (39) 
then reads 



(101) 



It constrains the frame {t^,z^) through V' in Eq. (87) 



To test the dynamical stability of the straight and 
static 'background' solution thus defined, we Fourier ex- 
pand the first order perturbations of the equations of 
motion, and derive the characteristic modes. If all the 
characteristic modes are real, the unperturbed solution 
is stable, otherwise it is unstable. 



A. Stability of the transverse modes 

As explained earlier, there are only two transverse de- 
grees of freedom and the corresponding equations are 
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given by Eq. ( 72 ) , which can be rewritten as 



in the general case of N currents. 



±/ ( u'^^^uP - ^v'^y.vP ] = 0, (102) 



where _L,/ is the projector orthogonal to the string world- 



sheet defined in Eq. (57 1. The string worldsheet is de- 



fined implicitly from its tangent space (generated by the 
eigenvectors u^^ and v^^). Eq. (61 1 therefore needs to be 
satisfied by the eigenvectors. It yields 



(103) 



Transverse perturbations of di^ do not change their 
mutual scalar products —Xij to first order and thus, the 
energy density U and tension T of the string are not 
affected by transverse perturbations of the worldsheet. 
Conversely, the current conservation equations are 
unaffected to first order by transverse perturbations. 
Thus, the transverse modes are decoupled and it is 
sufficient to consider transverse perturbations of the 
eigenvectors and u^, together with Eqs ( 102|103 l. 



This system of equations is exactly the same as in the 
case of a string carrying a single current [43j HH |45l [46] 



Given the usual definition [62] 5 (Q) = e''^e {Q) with 
— kfj^, k — k^yP, uj = —k^u^, where Q denotes 
uP or vP and e (Q) ^ 1, one finds that transverse per- 
turbations propagate at velocity 



T 



(104) 



as can also be checked directly from Eqs. p^02n and ( 103 1 



This result was to be anticipated, since the conservation 
equations are identical to those of the single string case 
and their projection orthogonal to the worldsheet do not 
depend on perturbations of the conserved currents but 
only on S (uP) and 5{vP). Thus, the dynamical stabil- 
ity of the transverse perturbation of the straight string 
requires simply that the tension be positive. 



T > 0. 



(105) 



In order to describe these modes, we once again 
assume that the string remains straight and the problem 
is effectively two-dimensional in the {tP,zP) plane. In 
this case, the unknowns are the string currents df. 
The corresponding equations can be chosen to be the 



conservation equations ( 69 1 and ( 75 1 , together with the 
equations of state (85 1. As in the single current case. 



the worldsheet is not perturbed by longitudinal modes. 
The tensors rjP''' and eP'^ are therefore not affected by 
longitudinal perturbations. 

We expand the perturbations of the N currents in 



Fourier modes as Sdif^e 



. The equations of motion 
for the longitudinal perturbations are then derived from 
Eq. (l69l) as 



r^iuj^eP^Sdj,) - 2- r^P-Sd^pdi^icj^eP^d^,) = 0, 

OXkl 



and from Eq. ( 75 1 as 



ifuj^Sd,, = 0. 



(106) 



(107) 



At this order in perturbations, a mode associated with 
one of the currents thus couples to all other unperturbed 
currents but not to their perturbated part. Thus, the 
string will be stable against longitudinal perturbations 
as long as the dispersion relation for a given mode is a 
real. Unfortunately, the system of 27V equations given 
by the two expressions above is quite complicated and 
there is a priori no way to solve it in the general case. In 
particular, as announced at the begining of this section, 
the simple form = —dT/dU > is not recovered. The 
main simplification which can be achieved is to halve the 
number of equations and unknowns by noting that, from 
Eq. (107), all the perturbations Jdj^ must be orthogonal 



to Ljfj,, and thus colinear to each other. However, this 
is not enough to solve the system algebraically or even 
to reduce the problem to a standard eigenvalue problem. 
Thus, to get explicit constraints of the elastic domain as- 
sociated with the longitudinal equations, we shall restrict 
ourselves to the case of two currents. 



B. Stability of the longitudinal modes 

We now consider the longitudinal modes, i.e., those 
propagating in an internal way along the string world- 
sheet. We expect the results to differ from those of 
the single string case since, contrary to the transverse 
conservation equations, the longitudinal conservation 
equations depend on perturbations of the conserved 
currents. The presence of more than one current implies 
current-current interactions. In addition, contrary to 
the single current case in which TP'^ is automatically 
diagonal, we loose the freedom to align the currents with 
our prefered frame {uP,vP). As a result, the well-known 
result = —dT/dU > is not expected to be recovered 



VII. THE ELASTIC DOMAIN FOR THE TWO 
CURRENTS MODEL 

We shall now concentrate explicitly on the elastic 
string model that pertains to the microscopic model of 
the first sections. In this case, the Lagrangian (64 1 that 



stems from Eq. ( 33 1 depends only on the state parame- 
ters Wi = Xii {i — 1,2, with the identification (p 1 and 
CT ^ 2), but not on the off-diagonal component x. Note 
that in the general case a dependence of C on xi2 = 2: is 
a priori allowed. In the case at hand, while Eq. ( [39| or 
Eq. (70) do not depend on a;, Eq. (87) does depend on 



it, and one must be careful to take this dependence into 
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account. The other relevant parameters of the problem, 
and /Cj, can similarly be identified with the ones de- 



fined in Eqs. (29 1 and (32) respectively with the change 
of indices (i — > 1 and a ^ 2. 



A. Stability of the transverse modes 

The propagation of transverse perturbations along the 
string takes place, independently of the string's internal 
structure, with a velocity il HSj SB] = T/U. 
Given that [/ > 0, the stability of a string subjected to 
transverse perturbations requires T > 0. This is equiva- 
lent to a constraint on x, which can be obtained in the 



microscopic theory using Eqs ( 46 1 , ( 47 1 and ( 48 1 , and 
reads 



X ^ ^lim 



1 



1 



-I 1/2 



ao8) 

This inequality is not automatically possible, for it de- 
mands > - (w'^IC^ + w^/C^) , a condition which has 
no counterpart in the single current case. 



B. Stability of the longitudinal modes 



The four longitudinal equations can be chosen as either 
the conservation of the four currents ^ and ^ i 
as the conservation of the two currents 2 ^^i^ of the 
longitudinal part of the stress energy tensor. Indeed, 
from Eq. ( [82| we know that the conservation of the four 
currents implies the conservation of the longitudinal part 
of the stress energy tensor, but also that the converse is 
true, since 

(109) 



Of course, because of the duality relations ( 90p2 1, we 
can symmetrically choose to keep the conservation of 
the two currents di 2^ and of the longitudinal part of the 
stress energy tensor instead. 



As in section VI B 



we work with the four conservation 



equations ( 69 1 and ( 75 ) together with the identity ( 86 ) 



and the equations of state ( 85 1 . The latter reduce to two 



equations only since the Lagrangian ( 33 1 only depends 
on the two state parameters wi^2- The corresponding 
dynamical equations perturbed to first order are then 
given by Eqs. ( |106|107[ ). 



1. General case 

We first will consider the cases when both fields are 
condensed in the string and a; 7^ or W1W2 7^ 0. The 
other cases, corresponding to having only one field 



condensed in the string and one field set to zero or 
having all four currents colinear and lightlike, will be 
studied separately. 

In the case at hand, the conserved currents have both 
timelike and spacelike components in the basis {u'^,v^) 
and it is therefore convenient to introduce the following 
two lightlike vectors 



1 



[d,^ - (-l)V/d..] 



(110) 



where the overall sign of e^'' is chosen such that — di^ 
if a; 7^ and 1011102 = 0. Another two lightlike vectors are 
chosen such that 



1 



0. 



2 ' ' if-t lu 

When W1W2 7^ 0, these vectors are given explicitly by 
1 

~ 2w. 



(Ill) 



(112) 



These four lightlike vectors are obviously related since 
there are only two null directions in the longitudinal 
plane. In fact, one can show that 



wid2p = xdi^, - sign(x) \/a;2 - wiW2{e^'' di^), (113) 
and 

1^2^!^ 



xd2i, -t- sign(a;) - wiW2(e^/d2i/), (114) 



where the sign of e^'^ was chosen appropriately when 
W1W2 7^ 0, and the sign function was used. Note that the 
term under the square root is always positive since, from 
Eq. (54 1, when W1W2 > 0, one must have x^ > W1W2 > 0. 



When a; = (and thus W1W2 < 0), the sign of x plays 
no role in these identities and can be fixed arbitrarily at 
-t-1. From Eqs ( |113 1 and (114) one can then derive the 
simple relation 



2fi 



with 



\ = X + sign(a;) \J x'^ — W1W2, 



(115) 



(116) 



Note again that A 7^ in the case we are considering. 
With these definitions, the currents can be rewritten as 



di. 



(117) 



Perturbations can be expanded in Fourier modes lo^ and 
are written as 



5d,,^{5dtel + 5dTe-\e^-^^ 



(118) 



with no summation on the repeated index i. The varia- 
tion of the state parameters is given by Eq. ( 84 1 as 



{5d^ +wMty 



(119) 
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The other two currents c*^ can be deduced from the iden- given by Eq. (1191. Note that since we assumed that 



tity (86) and 



two currents were condensed in the string, /C** ^ 0. 



^ = e^'^M,, + ^-^^Sw.e^'^d,,, (120) 



The equation of motion for the longitudinal modes are 
given by Eqs. (106' and (1071, and yield a homogeneous 



with no summation on the repeated index i, while bw^ linear system of equations. 



Jsr+<5d+ + Xf5dJ- = 0, 

+ X2"(5d2" = 0, 

{Xt + - wiXf + [-Xf + - wiXf )] Ml + ii.2(Xi+ - wxXl){w2Ml + (5d^) = 0, 

L2,l(^2^ - W-iX:^)(wxMX + + [^2^ + ^2,2(^2+ - W2X:,)W^ dd+ + [-X^ + L2.2{X+ - W2X^)\ Sd^ - 0, 

(121) 



where 



down in a compact form as 



^l-J — 



91n(/C") 



(122) 



The dispersion equation is obtained by equating the de- 
terminant to zero and is a homogeneous fourth degree 
polynomial in the two components of lo^ which does not 
appear to have any obvious solutions. It can be written 



Y,c,{x^f-\xiy = Q. 



(123) 



i=0 



Noting that XX2 = -^^2 1 from Eq. (1151, the coefficients 
Ci are found to be 



Co 

Cl 
C2 
C3 
C4 



= ~2\{Dwiwl + Dw2\^ + >? Li^i + wlL2,2) 

= D{wlwl + 4A^wiW2 + A'') + 4A^(1 + Lijwi + ^2,2^2) 



-2A (L'w?W2 + DwiX^ 



(124) 
(125) 
(126) 
(127) 
(128) 



where D is the determinant of the matrix Lij and A was 



defined in Eq. (116 1. We define 



X2 



-Qe. 



2t 



where C — ^t/^z is the dispersion relation. The string 
will remain stable under perturbations in the range of 
parameter space in which C is a real. Since the are 
real, it suffices that ^ be real. In terms of ^, Eq. (123) 
reads 



Co + ciC + C2r + Car + C4r = o. 



(130) 



Algebraic solutions to this equation are easily obtained 
using the well-known Ferrari procedure. Only and 
^2,2 J which depend implicitly on the w^'s, must be de- 
termined numerically from the integrated quantities of 



Section |IV[ In order to determine the regions of param- 
eter space in which the string is indeed stable, one must 
therefore input the values of Li^i and £2,2 obtained nu- 



(129) merically in the analytic forms of the roots of (130) ob- 



tained with the Ferrari procedure. We do so in the range 
w^min < < Winax ^ud —Xyi^ < X < xjim whcrc the Up- 
per and lower limits on wt are given by Eqs (|21| and ( 23 1 



and the limits on x are given by the condition for stabil- 
ity of the transverse modes. Figs [6] [7] [8] and [9] provide a 
comparison of the regions of stability obtained using the 
constraint x\\^ and the solutions of Eq. ( |130p G M for the 
cases of zero, weak and moderate coupling g. The region 
of stability in the a;-direction is given by the conditions 
for stability of the transverse modes while that in the 
Wi {i — 1,2) directions is provided by requiring that the 
roots of Eq. (130) be real. As shown in the figures, sta- 
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bility can be enhanced in the x direction as one goes to 
stronger couphng, but it is always significantly reduced 
in the Wi and W2 directions. This is especially true when 
one of the currents is spacelike: in Figs |7] and |8] the string 
is entirely unstable for a coupling greater or equal to ^3. 
Note that the region of stability of the analytic model 
of Section |VIII| is also shown in these figures. It is dis- 
cussed at the end of the next section. Finally, we point 
out once again that the stability conditions obtained by 



given by 



solving Eq. ( 130 1 is not valid when one or both currents 
are lightlike, i.e., for Wi'W2 = 0. The stability in these 
cases is discussed in the next subsections. 



2. Special case of a static condensate 

This case is characterized hy x — and either wi — 
or W2 = 0. It cannot be included in th e general case 



because A = and thus, from Eq. (1151, = 0. Let 
us call generically i the index for which Wi = 0, and j 
the other index. This case corresponds to uji = ki = 
in Eqs ^ and (10 1, and thus to a static condensate. As 



a consequence, it is clear that = 0. Using Eq. (106) 



one finds that the perturbation equations for the two cur- 
rents decouple. The two characteristic modes for the 
static condensate are simply given by the two light-like 
directions in the longitudinal plane. Thus, these pertur- 
bations are always stable. The stability condition for the 
other current reduces to 



1 



> 0. 



(131) 



It is interesting to note that the characteristic polyno- 
mial (123) found in the previous general case has two 



obvious real solutions corresponding to the two lightlike 
directions = and = 0, and two other solutions 



which are real exactly when the condition (131) is satis- 
fied. Thus, the general stability condition from the first 
case extends also to this limiting case. 



dT 
dU 



(l + 2u;,L,,,)-s"("''), (132) 



The stability condition is simply dT/dlJ < or equiva- 
Icntly that 1 -t- 2wiLi^i > 0. 



4- Special case of two static condensates 

This occurs when x — wi — W2 — 0. In this case, there 
is no current in the string, from (48 1, — = 0, thus 
from Eqs (45 1 and (46 1 U — T, and the stress energy 



tensor is Lorentz invariant. 

This limit corresponds effectively to a currentless 
Nambu-Goto string which is stable and although it is 
an interesting model in it own right, it is not an elastic 
string model. 



VIII. STABILITY OF THE SIMPLIFIED MODEL 

In this section, we test the possibility that while both 
fields are condensed in the string, they can be assumed 
(in some approximation) to be independent of each other. 
In this case the Lagrangian reads 



C = Ciixi) + C2{X2) 



(133) 



where Ci{Xi) is an effective Lagrangian which can be 
identified with that of a string with a single condensate. 
We include an additional term in order to compen- 
sate the bare Goto-Nambu energy contribution which, as 
we shall see, comes in twice, i.e., once for each of the 
Lagrangians £1 and £2- As shown in Refs HTJIH], the 
energy Ui and tension Ti as a function of the state pa- 
rameter Wi of a string with a single current are accurately 
reproduced by a Lagrangian of the following form 



Special case of only one current condensate 



In this case, the effective Lagrangian (33) depends 
only on one of the state parameters. There are only 
two longitudinal degrees of freedom and the equations 
derived in Subsection 1 of Section |VII B| do not ap- 
ply. Calling generically i the index of the remaining 
condensate, only Wi, K," and Li^i are nonzero. The 
perturbation equations derived from Eqs. ( |106 1 are 
exactly the same as the ones derived for the decoupled 
nonstatic condensate in the previous case. Thus the 



stability condition is also given by the inequality (131 1 



i.e., in the magnetic case, and 



ml 
2 



In 1 



Wi 



if Wi > 0, 

(134) 



if 



Wi < 0, 



(135) 

i.e., in the electric case. In Eqs. (134) and (135), m and 
rrii are adjustable mass parameters of the order of the 
Higgs mass nih = V^i] and of the order of the condensate 
mass parameters and ma- respectively. The energy 
per unit length Ui and tension Ti are given by [43l HH 



Note that this case has been studied in detail in previ- 
ous works 1331 1331 133 BE] and it is in fact well known that 
the speed of propagation of longitudinal perturbations is 



T, 



and 



if 



w^> Q 
(136) 
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and 



model. On the other hand, better agreement with the 
numerical results of the interacting model was found 
by keeping the effective scalar field mass parameters 
•nii equal to the best fit value obtained in the single 
condensate case. As expected, agreement of the model 
with the numerical results is very satisfactory for gi and 
92 but is unsatisfactory for 53. This serves to confirm 
that the model can only be used at weak coupling. We 
further note that (somewhat coincidentally) the model 
appears to better agree with the weakly coupled case 
(i.e., for ^2) than with the uncoupled case (i.e., for 
gi). Finally, we stress that the simplified model and 
the fully interacting model but with ^1 = should not 
be taken to be two descriptions of precisely the same 
physics since in the latter, there is an interplay between 
the two scalar fields through their couplings to the 
Higgs and associated C/(l)'°'^'*' gauge field. As a result, 
one should always expects a discrepency between the two. 

We now turn to a more detailed discussion of the stabil- 
ity of transverse and longitudinal modes in the simplified 
the conserved currents cf and depend only on the model. The vectors and are respectively timelike 



and T, = —C if Wi < 0. 

(137) 

A comparison between the numerical integration for a 
string with a single condensate and the analytic model 
is provided in Fig. [5] and shows remarkable agreement 
in most of the stable range given by Eq. (131) and the 
condition that T > 0. In practice, for the simplified 
model of Eq. (1331, i.e., in the situation where more 
than one field condensates onto the string, the mass 
parameters are shifted from the mass parameters one 
would get for a string with a single condensate. We 
will assume in this section that for the two fields to 
condensate, a minimum condition is that the associated 
currents be dynamically stable separately, i.e., ^ — ^ 
for each g31 HH HSJ SB]. Now, assuming that both 
fields condensate in the string, one can consider the dy- 
namical stability of this simplified string model against 
perturbations in the associated currents. Because of the 
decoupled form of the Lagrangian (1331, it is clear that 



^i{Xi)- This means that the longitudinal perturbation 
equations decouple and reduce to twice (once for each 
current) the usual longitudinal equations for a string 
with only one current. Thus, the stability condition is 
just cf^ > 0, a condition which we took to be necessary 
to have two current condensates. This result can be 
verified from Eqs. ( 121|122| by noting that for the 
Lagrangian o f Eq . (1331, Li_2 — ^2,1 
given by Eq. ([T32|. 



The stability of the string against transverse perturba- 
tions is simply given by (105 1. The stress energy tensor 
associated with Eq. (1331 is the sum of the stress energy 



tensors of each component 



and spacelike unit vectors associated with current i. In 
ord er to diagonalize T'"^, we write the stress energy ten- 
sor (138) in the basis {u^, ). The eigenvalues Ui and Ti 
are related by the equations of state in the usual way |37] 
of a one condensate string, and the other eigenvectors 
(^2, V 2) can be determined from (u^^ ) and x with Eqs 
(113) and (114). As stated before, the stability of the 



transverse modes is given by T > and this constraint 
on T can be used to obtain a constraint on the range 
of X. This is what we do in the remainder of this work. 
The longitudinal modes decouple and their stability, as in 
the single current case, is given by 1 -t- 2uJiLi^i > 0. For 



(138) 



the models of Eq. (134) and Eq. (135) this reduces to 
—7 < w < j/3 (or —m~ < Wi < mf/3) in both possible 
cases, i.e., wiW2 > and 'W1W2 < 0. 



where T/"^ are the stress energy tensors for a string with 
a single condensate. 



(139) 



in which Ui and Ti are the energy density and tension in 
the string. 

The energy density and tension of the simplified model 
for the two condensate case are shown in the lower part 
of Fig. [5] and are compared to the numerical results 
obtained for the interacting model in Figs. [3] and |4] for 
the three values of the coupling g considered in this work. 
In all three cases, the effective Higgs mass parameter m 
of the model (dashed lines) was adjusted so as to best 
fit the numerical results obtained for the full interacting 



5. Case I: W1W2 > 



In this case, > ^JwiW2, and = ^J~WiU^ for 
timelike currents and df^ — ^Jwlvi for spacelike currents. 



Using Eq. (113), one has 



(xu't + sign (x) \J x^ — w\W2v't \ , 



f sign (x) \/ x'^ — WiW2u'^ + xv'{\ 



(140) 



yJWlW2 



(141) 



Using Eq. (138), the stress energy tensor in the basis 



(m^, v'^) is given by the sum T^" -f m^rj'^'^ where 
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FIG. 5: (Top) Analytic models vs. numerically determined \J and T for a string with one current. The microscopic parameters 
used to compute I] and T numerically from the complete interacting case (solid lines) are those used for $ (left) and E (right) 
throughout the rest of this work. The mass parameters m and rrii (see Eqs. (134 1 and (135l) used to compute U and T in the 



analytic model (dashed lines) were adjusted to match the normalization of U and T obtained numerically. (Bottom) Analytic 
computation of U (left) and T (right) for a string with two currents using the mass parameters m as in the single current case 
shown in the upper part of the figure. 



Ui H 

W1W2 

WIW2 



{U2 - T2) 



\x\\/ x'^ — W1W2 \ 



V 



{U2 - T2) 



W1W2 



Ti 



X'T2 - (X-' - W1W2) U2 

WIW2 



(142) 



In order to address the first stability condition, T > 0), we compute U and T. They are simply given by the eigenvalues 
of the stress energy tensor. We therefore solve 



det [T'"' 



A) T?^''] ^X^ + [Ui + U2 + T1 + T2)X + 



{Ui + U2){Ti + T2)- 



V 



W1W2 



(Ui - Ti) {U2 - T2) 



= 0, 
(143) 



where X — X — m? and = x^ — w\W2 > 0. In order to 6. Case II: wiW2 < 

satisfy T > 0, one finds that 

, 2\frT^ rp 2\ ^^^^ case, and without loss of generality, we take 

y2 < y,^ y,, [Ui + U2-m )[I^ + l2-m ) ^ ajjd 4 = J^iJ2vli. We then have 

^ - ' ' (C/i-Ti)([/2-T2) ■ I i 1 ^ 



U9 = - fsign {x) \/ x'^ ~ WiW2U^ + xv'f \ , 



^/-WiW2 



21 



(145) 



■WlU'2 



^xw^ + sign (a;) 



(146) 



Following the procedure outlined in the preceding sec- 
tion, we find 



< \wiW2\ 



(C/i - Ti) ((72 - T2) 



(147) 



Alternatively, from Eq. ( 144[ ) and Eq. (1471, stable con- 
figurations can be expressed as a limit on x. For W1W2 > 
0, 



2 _wiW2{Ui-m^ + T2){U2-m^ + Ti) 

_ ^lim 



(C/i - Ti) (C/2 - T2) 



and for k;iW2 < 0, 



(148) 



X < X\u 



\wiW2\ {Ti + T2- m^) [Ui + U2- m^) 
(C/i-ri)(f/2-T2) 



(149) 

Figs. |6] [71 |8]and|9]provide a comparison of the regions of 
stability of the fully interacting model obtained using the 
constraint given by Eq. (108) and the solutions of 



Eq. ( 130 ) G M with the region of stability obtained for the 



simplified model using the constraints given by Eqs ( [148^ 
and ( 149 1 for gi, g2 and g^. In the simplified model, both 
Ti and Ui are proportional to m^. As a result, a;iiin is 
independent of m. It is interesting to note that, at least 
in principle, the coupling g indirectly plays a role in the 
physics of the decoupled model if the best fit m; 's differ in 
going from gi to ^3. If, however, the m^'s of the analytic 
model are kept the same for all values of g (which is the 
case here), the stability of the simplified model will be 
strictly the same independently of the value of g. In this 
case, the simplified model will reproduce the stability of 
the interacting full theory only for weak coupling. This is 
indeed the conclusion that can be drawn from the figures, 
where the regions of stability computed using simplified 
model approach those of the fully intereacting model. 



IX. CONCLUSIONS 

In the Witten bosonic neutral superconducting string 
model, the conserved current is given by the phase gradi- 
ent of a single [/(l)siobai s(2alar field condensate. At the 
field theory level and also in the integrated macroscopic 
worldsheet formalism |32l HH HSl US] , such a string can 
be described by a single state parameter given by the 
square of the phase gradient of the neutral field once the 
microscopic mass parameters and couplings are fixed. 
In this work, we extended both the existing microscopic 
and macroscopic descriptions of the neutral string to the 
case of a string carrying several of these currents. 



In the first part of this paper, we described how to 
generalize the single-current microscopic structure of a 
string to the two-current case (i.e., to the case of a string 
in the presence of two U{1)^^°^'^^ fields) and calculated 
the integrated quantities needed to describe the macro- 
scopic evolution of a string worldsheet endowed with two 
currents, i.e., the energy per unit length, the tension, and 
the total currents. We pointed out that at the field theory 
level, only the usual state parameters given by the square 
of the phase gradients are needed to fully describe the dy- 
namics. However, and in contrast with the other known 
case of a string carrying two currents |54j . in addition 
to the two Lorentz-invariant state parameters required 
in the microscopic description, the scalar product of the 
phase gradients of the different fields is needed in order 
to describe the energy per unit length, tension and in- 
tegrated currents relevant in the worldsheet description. 
This difference can be explained by noting that the field 
equations cannot depend on this third state parameter, 
denoted x, because x represents the (Lorentz) angle be- 
tween the worldsheet currents, and this could only have 
a dynamical effect on the microstructure if there existed 
an interaction between the scalar field phases. Such an 
interaction is excluded by the assumed symmetry. On 
the other hand, the macroscopic properties of the string 
needed to describe the worldsheet dynamics are non-local 
quantities and therefore can and do depend on the pa- 
rameter X. 

In this first part of the paper, we also solved nu- 
merically the field equations using a relaxation method 
and obtained the parameter range in which both fields 
condensate onto the string. As could be expected, 
when the C/(l)siobai f^gj^^jg strongly coupled at the 
microscopic level, the region of parameter space in which 
both fields condense onto the string is greatly reduced. 
In addition, we computed the energy per unit length, 
tension and the integrated current as a function of the 
first two state parameters. 

It was found that, at the microscopic level, the results 
obtained in the single current case extend to this more 
complicated case. In particular, the phase frequency 
threshold [H] that appears for a timelike current for 
one condensate exists here as well. As a result, in the 
region of the parameter space in which both currents 
play a significant part in the string physics, divergences 
are observed in the energy per unit length, tension 
and integrated currents. In particular, the divergences 
observed as one (but not the other) state parameter 
reaches its phase frequency threshold are similar to the 
ones observed in the single condensate string. However, 
a divergence of the tension T to positive values is 
observed when the phase frequency threshold of both 
currents is approached. This feature does not exist in 
the single condensate string: as the phase frequency 
threshold in the single condensate string is approached, 
the tension diverges negatively. This feature enlarges 
the range of validity of the no-spring conjecture |631 164j : 
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FIG. 6: Region of stability for Wa ~ — 7cr/2 and < (i.e., for timelike currents) for all three values of g and compared with 
the region of stability obtained with the analytic model. 




FIG. 7: Same as Fig. [6] but for Wcr ~ 7ct/2 and > (i.e., for spacelike currents). Note that in the case §3, the string is 
unstable in this range of parameters. 



it would affect the electromagnetically supported loop 
configuration |65j and their long-range gravitationnal 
properties [66]. Finally, exact relations describing the 
dependence of the various quantities of interest on 
the third state parameter were obtained. This was 
made possible only because, as already mentioned, the 
parameter x is absent at the microscopic level. 

In the second part of this paper, we presented a gen- 
eral extension of Carter's worldsheet formalism to de- 
scribe the worldsheet dynamics of a string endowed with 
N condensates. We derived the conservation equations 



associated with the 2N string currents and obtained a 
duality relation identical to the one that exists in the one 
condensate string. We then worked out the stability con- 
ditions of an iV-condensate string. At first order in per- 
turbations, the criterion for stability of transverse modes 
of a string carrying N currents is identical to that of a 
string carrying a single current: the propagation speed of 
perturbations must be positive. It reduces to the same 
requirement as in the single condensate string, namely 
that the string tension be positive. This conclusion can 
be reached by inspection, simply by substituting in the 
extrinsic equations of motion of the string the ansatz for 
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FIG. 8: Same as Fig. [6] but for Wo- ~ 7ct/2 and < (i.e., for one timelike and one spacelike current). As in the previous 
figure, the case of is unstable in this parameter range. 





FIG. 9: Same as Fig. [sjbut for Wc, ~ — 7ct/2 and > (i.e., roles reversed with respect to the previous figure). 



a perturbation. At first order, perturbations in the lon- 
gitudinal modes do not couple. The string will therefore 
be stable if the propagation speed of these modes is pos- 
itive. However, for a general A^-condensate string, this 
condition does not reduce to dT/dU < as in the single 
condensate string. Instead, one obtains a set of coupled 
equations for the perturbations of the A'^ currents (a cou- 
pling of the perturbation of a current with index i to all 
other unperturbed currents with index j ^ i). Further- 
more, in the worldsheet formalism, there exists an addi- 
tional set of state parameters, namely the scalar product 
of all pairs of distinct phase gradients that, in the N = 2 



condensate string, are naturally identified with the pa- 
rameter X of the microscopic field model. These new pa- 
rameters naturally complement the ones obtained when 
going from one to two currents. 

Given that there exists a priori no procedure to 
solve the set of longitudinal equations in the case of A^ 
currents, we carried out the numerical analysis in an 
application to a two-condensate string. In this case, the 
condition for stability of the transverse modes, namely 
T > 0, can be turned into a restriction on the range of 
possible variations condition on the third state param- 
eter. The perturbation equations in the longitudinal 



24 



direction reduces to a set of four coupled equations 
in which case the condition for stability turns into a 
constraint on the solutions of a 4"^ order polynomial 
equation. This condition is that its roots, which are 
functions of the dispersion relation of a given mode, be 
real. Algebraic solutions of this equation were obtained 
using the well-known Ferrari procedure. Combining the 
constraint on x with the constraints on the roots of the 
^th order polynomial equation completes the stability 
analysis of the two-condensate string. The regions of 
stability were shown in Figs [§ [T] [U and (9] 

In the third part of this work, we investigated 
whether the study of decoupled fields whose integrated 
Lagrangian can be reproduced using simple approximate 
analytic formulae is suitable to reproduce the physics in 
the two-coupled fields case. In this approximation, the 
Lagrangian is simply given by the sum of two single- 
current string Lagrangians. As expected, we found 
that even in this approximation, the energy per unit 
length and tension depend on both the single-current 
state parameters and x. Nevertheless, even when the 
nonlinear coupling constant is set to zero in the fully 
interacting Lagrangian, the simplified model should 
not be taken as a description of precisely the same 
physics because in the fully interacting Lagrangian, 
there is an interplay between the two scalar fields 
through their couplings to the Higgs and associated 

2^ -J local gauge field. As a result, one should always 
expect a discrepency between the two. In short, the full 
dynamical evolution of a two-current string, and thus 
presumably of an iV— current string, is not expected 
to be fully reproduced by such a model. However, the 



ability to give an approximate description of the physics 
in the interacting case with a fully analytic model is 
attractive and as shown in Figs |3] and |4] the energy per 
unit length and tension are very well approximated by 
the analytic description of uncoupled fields as long as 
the coupling g is small. Larger deviations are however 
observed for stronger couplings. 

We completed the study of the approximate model 
with an analysis of its stability against transverse and 
longitudinal perturbations in the currents. It was found 
that stability of the string against longitudinal pertur- 
bations yields the constrainst dT/dU < 0. This is of 
course identical to the result obtained for a string with a 
single condensate but different from the result obtained 
the second part of the paper in which coupled fields were 
considered. On the other hand, the stability of the string 
against transverse perturbations yields a constraint on 
X, as in the fully coupled case. We finally compared the 
regions of stability predicted by the analuytic model of 
uncoupled fields to the ones obtained in the interacting 
theory in Figs [6] |7] [8] and [9] and found good agreement 
at weak coupling only. We therefore conclude that it is 
possible to address the stability of string endowed with 
coupled fields in a fully analytic way with satisfactory 
accuracy at weak coupling only. 



Acknowledgments 

We thank B. Carter, G. Comer, A. C. Cordero and 
R. Prix for enlightening discussions. 



[1] T. W. B. Kibble. Topology of cosmic domains and 

strings. J. Phys. A, 9:1387-1398, 1976. 
[2] T. W. B. Kibble. Some implications of a cosmological 

phase transition. Phys. Rep., 67:183-199, 1980. 
[3] E. P. S. Shellard and A. Vilenkin. Cosmic strings and 

other topological defects. Cambridge University Press, 

Cambridge, England, 1994. 
[4] F. R. Bouchet, P. Peter, A. Riazuelo, and M. Sakellari- 

adou. Evidence against or for topological defects in the 

boomerang data? Phys. Rev. D , 65(2):021301, January 

2002. 

[5] P. Peter, F. R. Bouchet, A. Riazuelo, and M. SakeUar- 
iadou. Can topological defects still contribute to the 
cosmic microwave background radiation? In V. G. 
Gurzadyan, R. T. Jantzen, and R. Rufhni, editors. 
The Ninth Marcel Grossmann Meeting, pages 2013-2014, 
2002. 

[6] A. Linde. Axions in inflationary cosmology. Phys. Lett. 

B, 259:38-47, April 1991. 
[7] A. Linde. Hybrid inflation. Phys. Rev. D , 49:748-754, 

January 1994. 

[8] R. Jeannerot. Supersymmetric so(lO) model with infla- 
tion and cosmic strings. Phys. Rev. D , 53:5426-5436, 



May 1996. 

[9] R. Jeannerot. Inflation in supersymmetric unifled theo- 
ries. Phys. Rev. D , 56:6205-6216, November 1997. 

[10] R. Jeannerot, J. Rocher, and M. Sakellariadou. How 
generic is cosmic string formation in supersymmetric 
grand unified theories. Phys. Rev. D , 68(10):103514, 
November 2003. 

[11] G. Dvali, R. Kallosh, and A. Van Proeyen. D-term 
strings. Journal of High Energy Physics, 1:35, January 
2004. 

[12] S. C. Davis, A.-C. Davis, and M. Trodden. n=l super- 
symmetric cosmic strings. Phys. Lett. B, 405:257-264, 
February 1997. 

[13] T. Goto. Relativistic quantum mechanics of one- 
dimensional mechanical continuum and subsidiary con- 
dition of dual resonance model. Prog. Theor. Phys., 
46:1560-1569, November 1971. 

[14] Y. Nambu. Strings, monopoles, and gauge fields. Phys. 
Rev. D , 10:4262-4268, December 1974. 

[15] D. Austin, E. J. Copeland, and T. W. B. Kibble. Evo- 
lution of cosmic string configurations. Phys. Rev. D , 
48:5594-5627, December 1993. 

[16] A. Albrecht and N. Turok. Evolution of cosmic strings. 



25 



Phys. Rev. Lett, 54:1868-1871, April 1985. 

[17] A. Albrccht and N. Turok. Evolution of cosmic string 
networks. Phys. Rev. D , 40:973 1001, August 1989. 

[18] D. P. Bennett and F. R. Bouchet. Evidence for a scaling 
solution in cosmic-string evolution. Phys. Rev. Lett. , 
60:257-260, January 1988. 

[19] D. P. Bennett and F. R. Bouchet. Cosmic-string evolu- 
tion. Phys. Rev. Lett. , 63:2776-2779, December 1989. 

[20] D. P. Bennett and F. R. Bouchet. High-resolution sim- 
ulations of cosmic-string evolution, i. network evolution. 
Phys. Rev. D , 41:2408-2433, April 1990. 

[21] B. Allen and E. P. S. Shellard. Cosmic-string evolution: 
A numerical simulation. Phys. Rev. Lett. , 64:119-122, 
January 1990. 

[22] E. Witten. Superconducting strings. Nucl. Phys. B, 

249:557, 1985. 

[23] A.-C. Davis and S. C. Davis. Microphysics of so(lO) cos- 
mic strings. Phys. Rev. D , 55:1879-1895, February 1997. 

[24] R. L. Davis and E. P. S. Shellard. The physics of vor- 
tex superconductivity. Phys. Lett. B, 207:404-410, June 
1988. 

[25] R. L. Davis and E. P. S. Shellard. The physics of vortex 
superconductivity, ii. Phys. Lett. B, 209:485-490, August 
1988. 

[26] R. L. Davis and E. P. S. Shellard. Cosmic vortons. Nucl. 
Phys. B, 323:209-224, August 1989. 

[27] B. Carter. Cosmological relic distribution of conducting 
string loops. Ann. N. Y. Acad. Set., 647:758-766, De- 
cember 1991. 

[28] B. Carter. Mechanics of cosmic rings. Phys. Lett. B, 
238:166-171, April 1990. 

[29] B. Carter. How about vortons? In B. Guiderdoni, 
G. Greene, D. Hinds, and J. Tran Thanh van, editors. 
Dark Matter in Cosmology, Clocks and Tests of Funda- 
mental Laws, page 195, 1995. 

[30] R. Brandenberger, B. Carter, A.-C. Davis, and M. Trod- 
den. Cosmic vortons and particle physics constraints. 
Phys. Rev. D , 54:6059-6071, November 1996. 

[31] B. Carter and A.-C. Davis. Chiral vortons and cosmo- 
logical constraints on particle physics. Phys. Rev. D , 
61(12):123501, June 2000. 

[32] X. Siemens, V. Mandic, and J. Creighton. Gravitational- 
wave stochastic background from cosmic strings. Physical 
Review Letters, 98(11):111101, March 2007. 

[33] E. Komatsu ct al. Five-year wilkinson microwave 
anisotropy probe observations: Cosmological interpreta- 
tion. Ap. J. Supp., 180:330-376, February 2009. 

[34] N. Bevis, M. Hindmarsh, M. Kunz, and J. Urrestilla. Fit- 
ting CMB data with cosmic strings and inflation. ArXiv 
Astrophysics e-prints, February 2007. 

[35] A. L. Larsen. Dynamics of cosmic strings and springs - a 
covariant formulation. Class. Quantum Grav., 10:1541- 
1548, August 1993. 

[36] A. L. Larsen and M. Axenides. Runaway collapse of wit- 
ten vortex loops. Class. Quantum Grav., 14:443-453, 
February 1997. 

[37] B. Carter, P. Peter, and A. Gangui. Avoidance of collapse 
by circular current-carrying cosmic string loops. Phys. 
Rev. D , 55:4647-4662, April 1997. 

[38] A. Gangui, P. Peter, and C. Boehm. Could electromag- 
netic corrections solve the vorton excess problem? Phys. 
Rev. D , 57:2580-2589, February 1998. 

[39] B. Carter. Electromagnetic self interaction in strings. 
Phystcs Letters B, 404:246-252, February 1997. 



[40] P. Peter. Electromagnetically-supported cosmic string- 
loops. Physics Letters B, 298:60-62, Jarmary 1993. 

[41] X. Martin and P. Peter. Current-carrying string loop 
motion: Limits on the classical description and shocks. 
Phys. Rev. D , 61(4):043510, February 2000. 

[42] A. Cordero-Cid, X. Martin, and P. Peter. Current- 
carrying cosmic string loop 3d simulation: Towards a 
reduction of the vorton excess problem. Phys. Rev. D , 
65(8):083522, April 2002. 

[43] B. Carter. Covariant mechanics of simple and conduct- 
ing strings and membranes. In G. W. Gibbons, S. W. 
Hawking, and T. Vachaspati, editors, The Formation and 
Evolution of Cosmic Strings, page 143, 1990. 

[44] B. Carter. Mechanics of cosmic rings. Phys. Lett. B, 
238:166-171, April 1990. 

[45] B. Carter. Duality relation between charged elastic 
strings and superconducting cosmic strings. Phys. Lett. 
B, 224:61-66, June 1989. 

[46] B. Carter. Basic brane theory. Class. Quantum Grav., 
9:19, December 1992. 

[47] B. Carter and P. Peter. Supersonic string models for 
witten vortices. Phys. Rev. D , 52:R1744, August 1995. 

[48] A. Babul, T. Piran, and D. N. Spergel. Bosonic supercon- 
ducting cosmic strings, i. classical field theory solutions. 
Phys. Lett. B, 202:307-314, March 1988. 

[49] P. Peter. Superconducting cosmic string: Equation of 
state for spacelike and timelike current in the neutral 
limit. Phys. Rev. D , 45:1091-1102, February 1992. 

[50] P. Peter. Influence of the electric coupling strength in 
current-carrying cosmic strings. Phys. Rev. D , 46:3335- 
3349, October 1992. 

[51] Y. Lemperiere and E. P. S. Shellard. On the behaviour 
and stability of superconducting currents. Nucl. Phys. B, 
649:511-525, Jarmary 2003. 

[52] C. Ringeval. Equation of state of cosmic strings with 
fermionic current carriers. Phys. Rev. D , 63(6):063508, 
March 2001. 

[53] C. Ringeval. Fermionic massive modes along cosmic 
strings. Phys. Rev. D , 64(12) :123505, December 2001. 

[54] B. Carter. Cold, warm, and composite (cool) cosmic 
string models. Nucl. Phys. B, 412:345-371, January 1994. 

[55] L. D. Landau and E. M. Lifshitz. Course of Theoreti- 
cal Physics, Vol. 6: Fluid Mechanics, %130. Pergamon, 
Oxford, 1959. 

[56] P. G. Drazin and W. H. Reid. Hydrodyiianiic Stability. 
(Cambridge University Press, New- York, 1981. 

[57] N. Andersson, G. L. Comer, and R. Prix. Are pulsar 
glitches triggered by a superfluid two-stream instability? 
Phys. Rev. Lett. , 90(9):091101, March 2003. 

[58] T. Sidery, N. Andersson, and G. L. Comer. Waves and in- 
stabilities in dissipative rotating superfluid neutron stars. 
M. N. R. A. S., 385:335 348, March 2008. 

[59] D. Anderson, R. Fedele, and M. Lisak. A tutorial presen- 
tation of the two stream instability and landau damping. 
Am. J. Phys., 69:1262-1266, December 2001. 

[60] H. B. Nielsen and P. Olesen. Vortex-line models for dual 
strings. Nucl. Phys. B, 61:45-61, September 1973. 

[61] S. L. Adler and T. Piran. Relaxation methods for gauge 
field equilibrium equations. Rev. Mod. Phys., 56:1-40, 
January 1984. 

[62] B. Carter. Stability and characteristic propagation 
speeds in superconducting cosmic and other string mod- 
els. Phys. Lett. B, 228:466-470, September 1989. 

[63] P. Peter. No-cosmic-spring conjecture. Phys. Rev. D , 



47:3169-3176, April 1993. 

[64] P. Peter. Superconducting cosmic string tension: The 
positivity conjecture. In S. Gottlocber, J. P. Muecket, 
and V. Mueller, editors, Relativistic Astrophysics and 
Cosmology, page 320, 1992. 

[65] P. Peter. Electromagnetically-supported cosmic string- 
loops. Physics Letters B, 298:60-62, January 1993. 

[66] P. Peter and D. Puy. Gravitational field around spacelike 
and timelike current-carrying cosmic strings. Phys. Rev. 
D , 48:5546-5561, December 1993. 



